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MEASUREMENT OF GENE EXPRESSION PROFILES 
IN TOXICITY DETERMINATION 

5 Field of the Invention 

The invention relates generally to methods for detecting and monitoring 
phenotypic changes in in vitro and in vivo systems for assessing and/or determining 
the toxicity of chemical compounds, and more particularly, the invention relates to a 
method for detecting and monitoring changes in gene expression patterns in in vitro 
1 0 and in vivo systems for determining the toxicity of drug candidates. 

BACKGROUND 

The ability to rapidly and conveniently assess the toxicity of new compounds 
is extremely important. Thousands of new compounds are synthesized every year, 

1 5 and many are introduced to the environment through the development of new 

commercial products and processes, often with little knowledge of their short term 
and long term health effects. In the development of new drugs, the cost of assessing 
the safety and efficacy of candidate compounds is becoming astronomical: It is 
.estimated that the pharmaceutical industry spends an average of about 300 million 

20 dollars to bring a new pharmaceutical compound to market, e.g. Biotechnology, 13: 
226-228 (1995). A large fraction of these costs are due to the failure of candidate 
compounds in the later stages of the developmental process. That is. as the 
assessment of a candidate drug progresses from the identification of a compound as a 
drug candidate— for example, through relatively inexpensive binding assays or in vitro 

25 screening assays, to pharmacokinetic studies, to toxicity studies, to efficacy studies in 
model systems, to preliminary clinical studies, and so on, the costs of the associated 
tests and analyses increases tremendously. Consequently, it may cost several tens of 
millions of dollars to determine that a once promising candidate compound possesses 
a side effect or cross reactivity that renders it commercially infeasible to develop 

30 further. A great challenge of pharmaceutical development is to remove from further 
consideration as early as possible those compounds that are likely to fail in the later 
stages of drug testing. 

Drug development prograrrs are clearly structured with this objective in mind: 
however, rapidly escalating costs have created a need to develop even more stringent 

35 and less expensive screens in the early stages to identify false leads as soon as 

possible. Toxicity assessment is an area where such improvements may be made, for 
both drug development and for assessing the environmental, health, and safety effects 
of new compounds in general. 

- 1 - 

BNSDOOID:<WO . ., .9713877A1 I > 



3 



WO 97/13877 

PC77US96/16342 

mom.ormg the etTecs on . ^^"^ by 

* T blood chemisuy. weight ga L or T* 
5 .one. body temperature, re spi ra .ion rale letha i iIV ^1,-, T ^ * """^ 
provide a measure of the state of h Jl 7T ' """ Ch '""""'very 

such papers ^^0^ ""^ ^ *■« f 

compound Such "h T ^ ^ 3 ""^ ° f "» <*■ 

0 chemical on one occasion S^TS T " " adminiSM,i °" - 

be as long as a week or more In general „,„ Z „ """ S ° me » "»> 

-e .es, ehem.ca, on muhiple T <* 

or more „mes each day, iLuTaTas wh , ^ * one 

prolong .ypes of test, the •ta^S^^T^"^ ~ ^ 
- .s administered for a subsuntia, portion rf^^T " W, " Ch "* *" Chemical 
of .he mouse or ra, this is a period of 2 " 7™'' ~ 

7 years. years ' Ihe case °f the dog, it is for 5 to 

Significant costs are incurred in petokr u- 

■est animals for such ^ Z^ ZT & ^ ° f 

Moreover, because ofs^^ ^ ?" ^ 
ensure tha. a compouJTs free ^L" T ' ^ " "° ! 

however, provide a stajal H , hUmanS ' Suc " *»* *> 

of differentia.* efresTd "L^ „ T """" *** '° "» -I 
« I Science. 252: ,6^6* .991 m k " ^ ° r ~ «* Ad ™s 

of such apphcations indude increasing our knowledge of d^ale L ^"^ 
■denufying genes tha, p, ay importan , ro , es jn fc . 

d , e nos„c and therapeutic approaches tha, exp,o.. me ™ g^s"* 



WO 97/13877 



PCI7US96/16342 



products. While such approaches are attractive, those based on exhaustive, or even 
sampled, sequencing of expressed genes are still beset by the enormous effort 
required: It is estimated that 30-35 thousand different genes are expressed in a typical 
mammalian tissue in any given state, e.g. Ausubel et al, Editors, Current Protocols. 
5 5.8.1-5.8.4 (John Wiley & Sons, New York, 1992). Determining the sequences of 
even a small sample of that number of gene products is a major enterprise, requiring 
industrial-scale resources. Thus, the routine application of massive sequencing of 
expressed genes is still beyond current commercial technology. 

The availability of new assays for assessing the toxicity of compounds, such 
] 0 as candidate drugs, that would provide more comprehensive and precise information 
about the state of health of a test animal would be highly desirable. Such additional 
assays would preferably be less expensive, more rapid, and more convenient than 
current testing procedures, and would at the same time provide enough information to 
make early judgments regarding the safety of new compounds. 

15 

Summary of the Invention 
An object of the invention is to provide a new approach to toxicity assessment 
based on an examination of gene expression patterns, or profiles, in in vitro or in vivo 
. test systems. 

20 Another object of the invention is to provide a database on which to base 

decisions concerning the toxicological properties of chemicals, particularly drug 
candidates. 

A further object of the invention is to provide a method for analyzing gene 
expression patterns in selected tissues of test animals. 
25 A still further object of the invention is to provide a system for identifying 

genes which are differentially expressed in response to exposure to a test compound. 

Another object of the invention is to provide a rapid and reliable method for 
correlating gene expression with short term and long term toxicity in test animals. 

Another object of the invention is to identify genes whose expression is 
30 predictive of deleterious toxicity. 

The invention achieves these and other objects by providing a method for 
massively parallel signature sequencing of genes expressed in one or more selected 
tissues of an organism exposed to a test compound. An important feature of the 
invention is the application of novel DNA sorting and sequencing methodologies that 
35 permit the formation of gene expression profiles for selected tissues by determining 
the sequence of portions of many thousands of different polynucleotides in parallel. 
Such profiles may be compared with those from tissues of control organisms at single 
or multiple time points to identify expression patterns predictive of toxicity. 

-3- 
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the one or more tissues; and (g) correlating the frequency distribution of expressed 
genes in each of the one or more tissues with the toxicity of the compound. 

An important aspect of the invention is the identification of genes whose 
expression is predictive of the toxicity of a compound. Once such genes are 
5 identified, they may be employed in conventional assays, such as reverse transcriptase 
polymerase chain reaction (RT-PCR) assays for gene expression. 

Brief Description of the Drawings 
Figure 1 is a flow chart representation of an algorithm for generating 
1 0 minimally cross-hybridizing sets of oligonucleotides. 

Figure 2 diagrammatically illustrates an apparatus for carrying out 
polynucleotide sequencing in accordance with the invention. 

Definitions 

1 5 "Complement" or "tag complement" as used herein in reference to 

oligonucleotide tags refers to an oligonucleotide to which a oligonucleotide tag 
specifically hybridizes to form a perfectly matched duplex or triplex. In embodiments 
where specific hybridization results in a triplex, the oligonucleotide tag may be 
.selected to be either double stranded or single stranded. Thus, where triplexes are 

20 formed, the term "complement" is meant to encompass either a double stranded 

complement of a single stranded oligonucleotide tag or a single stranded complement 
of a double stranded oligonucleotide tag. 

The term "oligonucleotide" as used herein includes linear oligomers of natural 
or modified monomers or linkages, including deoxyribonucleosides, ribonucleosides, 

25 anomeric forms thereof, peptide nucleic acids (PNAs), and the like, capable of 
specifically binding to a target polynucleotide by way of a regular pattern of 
monomer-to-monomer interactions, such as Watson-Crick type of base pairing, base 
stacking, Hoogsteen or reverse Hoogsteen types of base pairing, or the like. Usually 
monomers are linked by phosphodiester bonds or analogs thereof to form 

30 oligonucleotides ranging in size from a few monomeric units, e.g. 3-4, to several tens 
of monomeric units. Whenever an oligonucleotide is represented by a sequence of 
letters, such as "ATGCCTG," it will be understood that the nucleotides are in 5'— >3' 
order from left to right and that "A" denotes deoxyadenosine, "C" denotes 
deoxycytidine, "G" denotes deoxyguanosine, and "T" denotes thymidine, unless 

35 otherwise noted. Analogs of phosphodiester linkages include phosphorothioate, 
phosphorodithioate, phosphoranilidate, phosphoramidate, and the like. Usually 
oligonucleotides of the invention comprise the four natural nucleotides; however, they 
may also comprise non-natural nucleotide analogs. It is clear to those skilled in the 
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As used herein, the term "complexity" in reference to a population of 
polynucleotides means the number of different species of molecule present in the 
population. 

As used herein, the terms "gene expression profile," and M gene expression 
5 pattern" which is used equivalently, means a frequency distribution of sequences of 
portions of cDNA molecules sampled from a population of tag-cDNA conjugates. 
Generally, the portions of sequence are sufficiently long to uniquely identify the 
cDNA from which the portion arose. Preferably, the total number of sequences 
determined is at least 1000; more preferably, the total number of sequences 
1 0 determined in a gene expression profile is at least ten thousand. 

As used herein, "test organism" means any in vitro or in vivo system which 
provides measureable responses to exposure to test compounds. Typically, test 
organisms may be mammalian cell cultures, particularly of specific tissues, such as 
hepatocytes. neurons, kidney cells, colony forming cells, or the like, or test organisms 
1 5 may be whole animals, such as rats, mice, hamsters, guinea pigs, dogs, cats, rabbits, 
pigs, monkeys, and the like. 

Detailed Description of the Invention 
The invention provides a method for determining the toxicity of a compound 

20 by analyzing changes in the gene expression profiles in selected tissues of test 
organisms exposed to the compound. The invention also provides a method of 
identifying toxicity markers consisting of individual genes or a group of genes that is 
expressed acutely and which is correlated with prolonged or chronic toxicity, or 
suggests that the compound will have an undesirable cross reactivity. Gene 

25 expression profiles are generated by sequencing portions of cDNA molecules 

construction from mRNA extracted from tissues of test organisms exposed to the 
compound being tested. As used herein, the term "tissue" is employed with its usual 
medical or biological meaning, except that in reference to an in vitro test system, such 
as a cell culture, it simply means a sample from the culture. Gene expression profiles 

30 derived from test organisms are compared to gene expression profiles derived from 

control organisms to determine the genes which are differentially expressed in the test 
organism because of exposure to the compound being tested. In both cases, the 
sequence information of the gene expression profiles is obtained by massively parallel 
signature sequencing of cDNAs, which is implemented in steps (c) through (f) of the 

35 above method. 

Toxicity Assessment 
Procedures for designing and conducting toxicity tests in in vitro and in vivo 
systems is well known, and is described in many texts on the subject, such as Loomis 
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apparatuses and chambers are available to perform studies for detecting effects of 
irritant or other toxic endpoints when they are administered by inhalation. The 
preferred method of administering an agent to animals is via the oral route, either by 
intubation or by incorporating the agent in the feed. 
5 Preferably, in designing a toxicity assessment, two or more species should be 

employed that handle the test compound as similarly to man as possible in terms of 
metabolism, absorption, excretion, tissue storage, and the like. Preferably, multiple 
doses or regimens at different concentrations should be employed to establish a dose- 
response relationship with respect to toxic effects. And preferably, the route of 

1 0 administration to the test animal should be the same as, or as similar as possible to, 
the route of administration of the compound to man. Effects obtained by one route of 
administration to test animals are not a priori applicable to effects by another route of 
administration to man. For example, food additives for man should be tested by 
admixture of the material in the diet of the test animals. 

1 5 Acute toxicity tests consist of administering a compound to test organisms on 

one occasion. The purpose of such test is to determine the symptomotology 
consequent to administration of the compound and to determine the degree of lethality 
of the compound. The initial procedure is to perform a series of range- finding doses 
jof the compound in a single species. This necessitates selection of a route of 

20 administration, preparation of the compound in a form suitable for administration by 
the selected route, and selection of an appropriate species. Preferably, initial acute 
toxicity studies are performed on either rats or mice because of their low cost, their 
availability, and the availability of abundant toxicologic reference data on these 
species. Prolonged toxicity tests consist of administering a compound to test 

25 organisms repeatedly, usually on a daily basis, over a period of 3 to 4 months. Two 
practical factors are encountered that place constraints on the design of such tests: 
First, the available routes of administration are limited because the route selected 
must be suitable for repeated administration without inducing harmful effects. And 
second, blood,- urine, and perhaps other samples, should be taken repeatedly without 

30 inducing significant harm to the test animals. Preferably, in the method of the 
invention the gene expression profiles are obtained in conjunction with the 
measurement of the traditional toxicologic parameters, such as listed in the table 
below: 
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assembled combinatorially from subunits which themselves are selected from a 
minimally cross-hybridizing set. For example, a set of minimally cross-hybridizing 
12-mers differing from one another by at least three nucleotides may be synthesized 
by assembling 3 subunits selected from a set of minimally cross-hybridizing 4-mers 
5 that each differ from one another by three nucleotides. Such an embodiment gives a 
maximally sized set of 9 3 , or 729, 12-mers. The number 9 is number of 
oligonucleotides listed by the computer program of Appendix la, which assumes, as 
with the 10-mers, that only 3 of the 4 different types of nucleotides are used. The set 
is described as ^maximal" because the computer programs of Appendices Ia-c provide 

1 0 the largest set for a given input (e.g. length, composition, difference in number of 
nucleotides between members). Additional minimally cross-hybridizing sets may be 
formed from subsets of such calculated sets. 

Oligonucleotide tags may be single stranded and be designed for specific 
hybridization to single stranded tag complements by duplex formation or for specific 

1 5 hybridization to double stranded tag complements by triplex formation. 

Oligonucleotide tags may also be double stranded and be designed for specific 
hybridization to single stranded tag complements by triplex formation. 

When synthesized combinatorially, an oligonucleotide tag preferably consists 
-of a plurality of subunits, each subunit consisting of an oligonucleotide of 3 to 9 

20 nucleotides in length wherein each subunit is selected from the same minimally cross- 
hybridizing set. In such embodiments, the number of oligonucleotide tags available 
depends on the number of subunits per tag and on the length of the subunits. The 
number is generally much less than the number of all possible sequences the length of 
the tag, which for a tag n nucleotides long would be 4 n . 

25 Complements of oligonucleotide tags attached to a solid phase support are 

used to sort polynucleotides from a mixture of polynucleotides each containing a tag. 
Complements of the oligonucleotide tags are synthesized on the surface of a solid 
phase support, such as a microscopic bead or a specific location on an array of 
synthesis locations on a single support, such that populations of identical sequences 

30 are produced in specific regions. That is, the surface of each support, in the case of a 
bead, or of each region, in the case of an array, is derivatized by only one type of 
complement which has a particular sequence. The population of such beads or regions 
contains a repertoire of complements with distinct sequences. As used herein in 
reference to oligonucleotide tags and tag complements, the term ki repertoire ?? means 

35 the set of minimally cross-hybridizing set of oligonucleotides that make up the tags in 
a particular embodiment or the corresponding set of tag complements. 

The polynucleotides to be sorted each have an oligonucleotide tag attached, 
such that different polynucleotides have different tags. As explained more fully 

- 11 - 
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1.30 x I0 9 


2.48 x I0 11 


62 


148 x I0 ? 


9.16x I0 8 


18 


105 x I0 5 


189 x JO 6 


39 


2.31 x !0 6 


9 02x I0 7 


332 


1.21 x I0 10 




28 


6.15 x JO 5 


172 x 10 7 


187 
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For some embodiments of the invention, where extremely large repertoires of 
tags are not required, oligonucleotide tags of a minimally cross-hybridizing set may 
be separately synthesized. Sets containing several hundred to several thousands, or 
5 even several tens of thousands, of oligonucleotides may be synthesized directly by a 
variety of parallel synthesis approaches, e.g. as disclosed in Frank et al, U.S. patent 
4,689.405; Frank et al, Nucleic Acids Research, 1 1 : 4365-4377 (1983); Matson et al. 
Anal. Biochem.. 224: 1 10-1 16 (1995); Fodor et al, International application 
PCT/US93/04 145; Pease et al, Proc. Natl. Acad. Sci., 91 : 5022-5026 (1994); 

10 Southern et al, J. Biotechnology, 35: 217-227 (1994), Brennan, International 

application PCT/US94/05896; Lashkari et al, Proc. Natl. Acad. Sci., 92: 7912-7915 
(1995); or the like. 

Preferably, oligonucleotide tags of the invention are synthesized 
combinatorial ly out of subunits between three and six nucleotides in length and 

1 5 selected from the same minimally cross-hybridizing set. For oligonucletides in this 
range, the members of such sets may be enumerated by computer programs based on 
the algorithm of Fig. I . 

The algorithm of Fig. 1 is implemented by first defining the characteristics of 
the subunits of the minimally cross-hybridizing set, i.e. length, number of base 

20 differences between members, and composition, e.g. do they consist of two, three, or 
four kinds of bases. A table M n , n=l, is generated (100) that consists of all possible 
sequences of a given length and composition. An initial subunit S \ is selected and 
compared (120) with successive subunits Sj for i=n+l to the end of the table. 
Whenever a successive subunit has the required number of mismatches to be a 

25 member of the minimally cross-hybridizing set, it is saved in a new table M n +] (125), 
that also contains subunits previously selected in prior passes through step 120. For 
example, in the first set of comparisons, M2 will contain Sj; in the second set of 
comparisons, M3 will contain Sj and S2; in the third set of comparisons, M4 will 
contain Sj, S2, and S3; and so on. Similarly, comparisons in table Mj will be 

30 between Sj and all successive subunits in Mj. Note that each successive table M n +j 
is smaller than its predecessors as subunits are eliminated in successive passes 
through step 130. After every subunit of table M n has been compared (140) the old 
table is replaced by the new table M n + j , and the next round of comparisons are 
begun. The process stops (1 60) when a table M n is reached that contains no 

35 successive subunits to compare to the selected subunit Sj, i.e. M n =M n +i . 

Preferably, minimally cross-hybridizing sets comprise subunits that make 
approximately equivalent contributions to duplex stability as every other subunit in 
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duplex stabilities, e.g. Rychlik et al Nu,!,- f P ™ WrS ^""i 

and I8: „o«4 12 <U ^^H^*'* 0989> 

W — • C* Biochem. Mol. Z ,, 2 T^^l 

For shorter tags. e.g. about 30 nucleotides or less the 7 I * "* ''^ 

and Wetmur is preferred, and for longer £ * 

an algorithm disclosed by Suggs e , a, pagcTo^ 6 „T « " " ^ 
Svmp. Dev. Bio,., Vol. 23 (Academic'^ ^ «*»• 

designing sets of minimalJycross-hvbridi 7 in 0 u • for 
■ermina. nucleotides when subunits are as^IL f ■«««"<*■»• energies of 




25 



30 



35 



where the primed W's indicate complements With a r 

Perfectly matched duplexes, all JLZTZ T^ ltZT^ 'T* 

thereby reducing the stabilitv «r , , *' m,sma *hes 

mismatched wo* att^^ ^ ^,7^ ^ ^ 

mismatches a. signiScan.lv less 1 
terminus. uupiexes with the same mismatch at a 

s^J^ZZ^Z* T m " y ™ - - those whose 

more ful.y below Z aL^Te f nUC ' e ° tidCS ' M d " 

exonudease activity of a DNA polvme^e 7 PP ° m * ° f "* 
minimally cross-hybridizing set of J»T T " S * " " m "'^ 

from the group JL^O^T *" ^ 
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Table II 

Word: w, w 2 w 3 w 4 

Sequence: GATT TGAT TAGA TTTG 

Word: w 5 w 6 w 7 w 8 

Sequence: GTAA AGTA ATGT AAAG 



In this set, each member would form a duplex having three mismatched bases with 
1 0 the complement of every other member. 

Further exemplary minimally cross-hybridizing sets are listed below in Table 
III. Clearly, additional sets can be generated by substituting different groups of 
'nucleotides, or by using subsets of known minimally cross-hybridizing sets. 



15 

Table III 

Exemplary Minimally Cross-Hybridizing Sets of 4-mer Subunits 



Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


Set 6 


CATT 


ACCC 


AAAC 


AAAG 


AACA 


AACG 


CTAA 


AGGG 


ACCA 


ACCA 


ACAC 


ACAA 


TCAT 


CACG 


AGGG 


AGGC 


AGGG 


AGGC 


ACTA 


CCGA 


CACG 


CACC 


CAAG 


CAAC 


TACA 


CGAC 


CCGC 


CCGG 


CCGC 


CCGG 


TTTC 


GAGC 


CGAA 


CGAA 


CGCA 


CGCA 


ATCT 


GCAG 


GAGA 


GAGA 


GAGA 


GAGA 


AAAC 


GGCA 


GCAG 


GCAC 


GCCG 


GCCC 




AAAA 


GGCC 


GGCG 


GGAC 


GGAG 
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Set / 

AAGA 

ACAC 

AGCG 

CAAG 

CCCA 

CGGC 

GACC 

GCGG 

GGAA 



Set 8 

AAGC 

ACAA 

AGCG 

CAAG 

CCCC 

CGGA 

GACA 

GCGG 

GGAC 



Set 9 

AAGG 

ACAA 

AGCC 

CAAC 

CCCG 

CGGA 

GACA 

GCGC 

GGAG 



Set 10 
ACAG 
AACA 
AGGC 
CAAC 
CCGA 
CGCG 
GAGG 
GCCC 
GGAA 



Set 11 
ACCG 
AAAA 
AGGC 
CACC 
CCGA 
CGAG 
GAGG 
GCAC 
GGCA 
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Set 12 
ACGA 
AAAC 
AGCG 
CACA 
CCAG 
CGGC 
GAGG 
GCCC 
GGAA 
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mixtures of the activated monomers to the growing oligonucleotide during the 
coupling steps. Preferably, oligonucleotide tags and tag complements are synthesized 
on a DNA synthesizer having a number of synthesis chambers which is greater than or 
equal to the number of different kinds of words used in the construction of the tags. 

5 That is, preferably there is a synthesis chamber corresponding to each type of word. 
In this embodiment, words are added nucleotide-by-nucleotide, such that if a word 
consists of five nucleotides there are five monomer couplings in each synthesis 
chamber. After a word is completely synthesized, the synthesis supports are removed 
from the chambers, mixed, and redistributed back to the chambers for the next cycle 

1 0 of word addition. This latter embodiment takes advantage of the high coupling yields 
of monomer addition, e.g. in phosphoramidite chemistries. 

Double stranded forms of tags may be made by separately synthesizing the 
complementary strands followed by mixing under conditions that permit duplex 
formation. Alternatively, double stranded tags may be formed by first synthesizing a 

1 5 single stranded repertoire linked to a known oligonucleotide sequence that serves as a 
primer binding site. The second strand is then synthesized by combining the single 
stranded repertoire with a primer and extending with a polymerase. This latter 
approach is described in Oliphant et ah Gene, 44: 1 77- 1 83 ( 1 986). Such duplex tags 
. may then be inserted into cloning vectors along with target polynucleotides for sorting 

20 and manipulation of the target polynucleotide in accordance with the invention. 

When tag complements are employed that are made up of nucleotides that 
have enhanced binding characteristics, such as PNAs or oligonucleotide NS'-^PS* 
phosphoramidates. sorting can be implemented through the formation of D-loops 
between tags comprising natural nucleotides and their PNA or phosphoramidate 

25 complements, as an alternative to the "stripping" reaction employing the 3'— >5" 
exonuclease activity of a DNA polymerase to render a tag single stranded. 

Oligonucleotide tags of the invention may range in length from 12 to 60 
nucleotides or basepairs. Preferably, oligonucleotide tags range in length from 1 8 to 
40 nucleotides or basepairs. More preferably, oligonucleotide tags range in length 

30 from 25 to 40 nucleotides or basepairs. In terms of preferred and more preferred 
numbers of subunits, these ranges may be expressed as follows: 



Table IV 

Numbers of Subunits in Tags in Preferred Embodiments 

35 

Monomers 

in Subunit Nucleotides in Oligonucleotide Tag 

(12-60) (18-40) (25-40) 
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3 «0«*»i«. 6-13 subunits 8-13 subunits 

4-10subunits 6-10subunits 
3-8 subunits 5-8 subunits 

3-6subunits 4-6 subunits 



4 3-15 subunits 

5 2-12 subunits 

6 2- 10 subunits 



i- '000 members; ^S^' T" 0 "- ° f ^ "•• — - a, 

1 0,000 m ernbers. '""^"'"•^"^"f^^conuina,,^, 

Triplex Tag c 
Genera.,* U,ird slraild associatio „ vja 

homopyrimidine-nomopurine Backs in a dTT ^ ° fb,nd "'« ls m °* ".We along 

Wica.es „„o g s,ee„ rypc B mJ n ^Z^TlZ iT'^T **** 
examp,e. HoogsIeen _ ^ ^ ^ ^XX^eT'^' 

and the like) in order ,o ^L Z !^ ™'™ ( °* """l""* C »°™- 

particular embodies 7 B 7oL ,7" re8U ' a,e ' """" SUbi,i ^ - **«> « 

Acad. Sci.. 93: 4320-4325 H996V n r X " Proc ' NaU - 

Chen, See., , , 4: f 4 ^ ^Beli " ^ « * J A ™ 

Giovanni e, a,, Proc. uJaZ Tv Z^ ' '* W6 " 4 ' 82 " 
Science, 238. 645.650(1987)- M Z' , ( '" 2); MoSerand De "an. 
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32: 666-690 (1993); Escude et al, Proc. Natl. Acad. Sci., 93: 4365-4369 (1996); and 
the like. Conditions for annealing single-stranded or duplex tags to their single- 
stranded or duplex complements are well known, e.g. Ji et al, Anal. Chem. 65: 1323- 
1328 (1993); Cantor et al, U.S. patent 5,482,836; and the like. Use of triplex tags has 
5 the advantage of not requiring a ^stripping" reaction with polymerase to expose the 
tag for annealing to its complement. 

Preferably, oligonucleotide tags of the invention employing triplex 
hybridization are double stranded DNA and the corresponding tag complements are 
single stranded. More preferably, 5-methylcytosine is used in place of cytosine in the 

1 0 tag complements in order to broaden the range of pH stability of the triplex formed 
between a tag and its complement. Preferred conditions for forming triplexes are 
fully disclosed in the above references. Briefly, hybridization takes place in 
concentrated salt solution, e.g. 1.0 M NaCI, 1.0 M potassium acetate, or the like, at 
pH below 5.5 ( or 6.5 if 5-methylcytosine is employed). Hybridization temperature 

1 5 depends on the length and composition of the tag; however, for an 1 8-20-mer tag of 
longer, hybridization at room temperature is adequate. Washes may be conducted 
with less concentrated salt solutions, e.g. 10 mM sodium acetate, 100 mM MgCl 2 , pH 
5.8, at room temperature. Tags may be eluted from their tag complements by 
- incubation in a similar salt solution at pH 9.0. 

20 Minimally cross-hybridizing sets of oligonucleotide tags that form triplexes 

may be generated by the computer program of Appendix Ic, or similar programs. An 
exemplary set of double stranded 8-mer words are listed below in capital letters with 
the corresponding complements in small letters. Each such word differs from each of 
the other words in the set by three base pairs. 

25 

Table V 

Exemplary Minimally Cross-Hybridizing 
Set of DoubleStranded 8-mer Tags 



Z. f 


-AAGGAGAG 


5' 


-AAAGGGGA 


5' 


-AGAGAAGA 


5' 


-AGGGGGGG 


3' 


-TTCCTCTC 


3' 


-TTTCCCCT 


3' 




3 ' 


-TCCCCCCC 




-ttcctctc 


3' 


-t ttcccct 


3' 


- tctctt ct 


3' 


- tccccccc 


^ * 


-AAAAAAAA 


S' 


-AAGAGAGA 


C r 


-AGGAAAAG 


5' 


-GAAAGGAG 


2 * 


-TTTTTTTT 


3' 


-TTCTCTCT 


3' 


-TCCTTTTC 


3' 


-CTTTCCTC 


3 ' 


-tltttttt 


3' 


-ttctctct 


3' 


-tccttttc 


3' 


-ctttcctc 


C ' 


-AAAAAGGG 


5' 


-AGAAGAGG 


5' 


-AGGAAGGA 


5' 


-GAAGAAGG 


'1 ' 


-TTTTTCCC 


3' 


-TCTTCTCC 


3' 




3' 


-CTTCTTCC 


1 r 
~> 


-tttttccc 


3' 


-tcttctcc 


3' 


-tccrtcc: 


3 ' 


-cttct:c: 




-AAAGGAAG 


b' 


-AGAAGGAA 


5' 


-AGGGGAAA 


5 ' 


-GAAGAGAA 


7 * 




3' 


-TCTTCCTT 


3' 




3 ' 


— CTTCTCTT 


3' 


-::tccttc 


3' 


-tcttcctt 


3' 


- 1 ccccttt 


3' 


-cttctctt 
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Table VI 



S^STt^ ° f VariOUS D ° Ub,e Stranded Tags 



Nucleotide 
Difference 

Oiio^f 6 " -, Maximal Size 

Oligonucieond « M ™™"> 



e 



of Minimally ".".nuu.y Sjzeof 

Cross- u..t rC l SS " Repertoire 



HvbriH c Hybridizing C Size of 

U aBL ^_«*"^S« Se , R 

rive Words 



4 

6 



8 
8 



4096 3.2 x,0 4 

8 3 4096 3.2 x,0< 

io , 16 65 *J° 4 1 .05 x I0 6 



i5 
20 
20 
20 



5 
5 
6 
8 
JO 



8 

92 
765 
92 
22 
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wide variety of compositions, including glass, plastic, silicon, alkanethiolate- 
derivatized gold, cellulose, low cross-linked and high cross-linked polystyrene, silica 
gel, polyamide. and the like. Preferably, either a population of discrete particles are 
employed such that each has a uniform coating, or population, of complementary 
5 sequences of the same tag (and no other), or a single or a few supports are employed 
with spatially discrete regions each containing a uniform coating, or population, of 
complementary sequences to the same tag (and no other). In the latter embodiment, 
the area of the regions may vary according to particular applications; usually, the 
regions range in area from several urn 2 , e.g. 3-5, to several hundred vm*. e.g. 100- 
1 0 500. Preferably, such regions are spatially discrete so that signals generated by 

events, e.g. fluorescent emissions, at adjacent regions can be resolved by the detection 
system being employed. In some applications, it may be desirable to have regions 
with uniform coatings of more than one tag complement, e.g. for simultaneous 
sequence analysis, or for bringing separately tagged molecules into close proximity. 
1 5 Tag complements may be used with the solid phase support that they are 

synthesized on. or they may be separately synthesized and attached to a solid phase 
support for use, e.g. as disclosed by Lund et al, Nucleic Acids Research, 1 6: 10861- 
10880 (1988); Albretsen et al. Anal. Biochem., 189: 40-50 (1990); Wolf et al. Nucleic 
- Acids Research, 15: 291 1-2926 (1987); or Ghosh et al. Nucleic Acids Research, 15. 
20 5353-5372 ( 1 987). Preferably, tag complements are synthesized on and used with the 
same solid phase support, which may comprise a variety of forms and include a 
variety of linking moieties. Such supports may comprise microparticles or arrays, or 
matrices, of regions where uniform populations of tag complements are synthesized. 
A wide variety of microparticle supports may be used with the invention, including 
25 microparticles made of controlled pore glass (CPG), highly cross-linked polystyrene, 
acrylic copolymers, cellulose, nylon, dextran, latex, polyacrolein, and the like, 
disclosed in the following exemplary references: Meth. Enzymol., Section A, pages 
11 -147, vol. 44 (Academic Press, New York, 1976); U.S. patents 4,678.814: 
4.413,070; and 4,046;720; and Pon. Chapter 19, in Agrawal, editor. Methods in 
30 Molecular Biology, Vol. 20, (Humana Press, Totowa, NJ, 1 993). Microparticle 
supports further include commercially available nucleoside-derivatized CPG and 
polystyrene beads (e.g. available from Applied Biosystems. Foster City, CA); 
derivatized magnetic beads; polystyrene grafted with polyethylene glycol (e.g.. 
TentaGelTM, Rapp Polymere, Tubingen Germany), and the like. Selection of the 
5 support characteristics, such as material, porosity, size, shape, and the like, and the 
type of linking moiety employed depends on the conditions under which the tags are 
used. For example, in applications involving successive processing with enzymes, 
supports and linkers that minimize steric hindrance of the enzymes and that facilitate 
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numbers of microparticles on planar supports, such as a microscope slide. 
Particularly preferred non-porous microparticles are the glycidal methacrylate (GMA) 
beads available from Bangs Laboratories (Carmel, IN). Such microparticles are 
useful in a variety of sizes and derivatized with a variety of linkage groups for 
5 synthesizing tags or tag complements. Preferably, for massively parallel 

manipulations of tagged microparticles, 5 \xxn diameter GMA beads are employed. 



10 

Attaching Tags to Polynucleotides 
For Sorting onto Solid Phase Supports 
An important aspect of the invention is the sorting and attachment of a 
populations of polynucleotides, e.g. from a cDNA library, to microparticles or to 

1 5 separate regions on a solid phase support such that each microparticle or region has 
substantially only one kind of polynucleotide attached. This objective is 
accomplished by insuring that substantially all different polynucleotides have 
different tags attached. This condition, in turn, is brought about by taking a sample of 
. the full ensemble of tag-polynucleotide conjugates for analysis. (It is acceptable that 

20 identical polynucleotides have different tags, as it merely results in the same 

polynucleotide being operated on or analyzed twice in two different locations.) Such 
sampling can be carried out either overtly—for example, by taking a small volume 
from a larger mixture—after the tags have been attached to the polynucleotides, it can 
be carried out inherently as a secondary effect of the techniques used to process the 

25 polynucleotides and tags, or sampling can be carried out both overtly and as an 
inherent pan of processing steps. 

Preferably, in constructing a cDNA library where substantially all different 
cDN As have different lags, a tag repertoire is employed whose complexity, or number 
of distinct tags, greatly exceeds the total number of mRNAs extracted from a cell or 

30 tissue sample. Preferably, the complexity of the tag repertoire is at least 10 times that 
of the polynucleotide population; and more preferably, the complexity of the tag 
repertoire is at least 100 times that of the polynucleotide population. Below, a 
protocol is disclosed for cDNA library construction using a primer mixture that 
contains a full repertoire of exemplary 9-word tags. Such a mixture of tag-containing 

35 primers has a complexity of 8 9 , or about 1 .34 x 10 8 . As indicated by Winslow et al, 
Nucleic Acids Research, 19: 3251-3253 (1991), mRNA for library construction can 
be extracted from as few as 10-100 mammalian cells. Since a single mammalian cell 
contains about 5 x 10 5 copies of mRNA molecules of about 3.4 x 10 4 different kinds. 
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10 8 ). then X=.00746 and P(2)=2.76 x 10' 5 . Thus, a sample of one million molecules 
gives rise to an expected number of doubles well within the preferred range. Such a 
sample is readily obtained as follows: Assume that the 5 x 10 1 1 mRNAs are perfectly 
converted into 5 x 10 n vectors with tag-cDNA conjugates as inserts and that the 5 x 
5 1 0 1 1 vectors are in a reaction solution having a volume of 100 |aK Four 1 0-fold serial 
dilutions may be carried out by transferring 10 |il from the original solution into a 
vessel containing 90 jal of an appropriate buffer, such as TE. This process may be 
repeated for three additional dilutions to obtain a 100 jal solution containing 5 x 10 5 
vector molecules per |il. A 2 |il aliquot from this solution yields 10 6 vectors 

1 0 containing tag-cDNA conjugates as inserts. This sample is then amplified by straight 
forward transformation of a competent host cell followed by culturing. 

Of course, as mentioned above, no step in the above process proceeds with 
perfect efficiency. In particular, when vectors are employed to amplify a sample of 
tag-polynucleotide conjugates, the step of transforming a host is very inefficient. 

1 5 Usually, no more than 1% of the vectors are taken up by the host and replicated. 

Thus, for such a method of amplification, even fewer dilutions would be required to 
obtain a sample of 1 0 6 conjugates. 

A repertoire of oligonucleotide tags can be conjugated to a population of 
- polynucleotides in a number of ways, including direct enzymatic ligation, 

20 amplification, e.g. via PCR, using primers containing the tag sequences, and the like. 
The initial ligating step produces a very large population of tag-polynucleotide 
conjugates such that a single tag is generally attached to many different 
polynucleotides. However, as noted above, by taking a sufficiently small sample of 
the conjugates, the probability of obtaining "doubles," i.e. the same tag on two 

25 different polynucleotides, can be made negligible. Generally, the larger the sample 
the greater the probability of obtaining a double. Thus, a design trade-off exists 
between selecting a large sample of tag-polynucleotide conjugates— which, for 
example, ensures adequate coverage of a target polynucleotide in a shotgun 
sequencing operation or adequate representation of a rapidly changing mRNA pool. 

30 and selecting a small sample which ensures that a minimal number of doubles wall be 
present. In most embodiments, the presence of doubles merely adds an additional 
source of noise or, in the case of sequencing, a minor complication in scanning and 
signal processing, as microparticles giving multiple fluorescent signals can simply be 
ignored. 

35 As used herein, the term "substantially all" in reference to attaching tags to 

molecules, especially polynucleotides, is meant to reflect the statistical nature of the 
sampling procedure employed to obtain a population of tag-molecule conjugates 
essentially free of doubles. The meaning of substantially all in terms of actual 
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percentages of tag-molecule conjugates depends on how the tags are being employed 
Preferably, for nucleic acid sequencing, substantially all means that at least eighty 
percent of the polynucleotides have unique tags attached. More preferably, it means 
that at least ninety percent of the polynucleotides have unique tags attached Still 
more preferably, it means that at least ninety-five pe rcent of the polynucleotides have 
unique tags attached. And, most preferably, it means that at least ninety-nine percent 
of the polynucleotides have unique tags attached. 

Preferably, when the population of polynucleotides consists of messenger 
(mRNA) ' ^nucleotides tags may be attached by reverse transcribing the 
mRNA w,th a set of primers preferably containing complements of tag sequences 
An exemplary set of such primers could have the following sequence (SEQ ID NO: 

5 ' -mRNA- [A] n -3' 

[T] 19 GG [W, W, W, C] qAC CAGCTG ATC- 5 • -biotin 



where "[W,W,W,C] 9 » represents the sequence of an oligonucleotide tag of nine 
- subunits of four nucleotides each and »[W,W,W,C]" represents the subuni, sequences 
listed above, i.e. »W» represents T or A. The underlined sequences identify an 
optional restriction endonuclease site that can be used to release the polynucleotide 
from attachment to a solid phase support via the biotin, if one is employed For the 
above pnmer, the complement attached to a microparticle could have the form: 

5 ' - [G, W, W, W] gTGG-linker-microparticle 

After reverse transcription, the mRNA is removed, e.g. by RNase H digestion 
and the second strand of the cDNA is synthesized using, for example, a primer of the" 
following form (SEQ ID NO: 2): 

5 1 -NRRGATCYNNN-3 • 

where N is any one of A, T, G, or C; R is a purine-containing nucleotide, and Y is a 
pynmidme-containing nucleotide. This particular primer creates a Bst Y 1 restriction 
site in the resulting double stranded DNA which, together with the Sal I site 
facilitates cloning into a vector with, for example. Bam HI and Xho 1 sites After Bst 
Yl and Sal I digestion, the exemplary conjugate would have the form- 
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5 1 -RCGACCA [C, W, W, W] 9GG [T] 29- cDNA -NNNR 

GGT[G,W,W,W] 9 CC[A] 19 - r DNA -NNNYCTAG-5 ' 

The polynucleotide-tag conjugates may then be manipulated using standard molecular 
5 biology techniques. For example, the above conjugate- which is actually a mixture- 
may be inserted into commercially available cloning vectors, e.g. Stratagene Cloning 
System (La Jolla, CA); transfected into a host, such as a commercially available host 
bacteria; which is then cultured to increase the number of conjugates. The cloning 
vectors may then be isolated using standard techniques, e.g. Sambrook et al, 
10 Molecular Cloning, Second Edition (Cold Spring Harbor Laboratory, New York ? 

1989). Alternatively, appropriate adaptors and primers may be employed so that the 
conjugate population can be increased by PCR. 

Preferably, when the ligase-based method of sequencing is employed, the Bst 
Yl and Sal I digested fragments are cloned into a Bam HI-/Xho I-digested vector 
1 5 having the following single-copy restriction sites (SEQ ID NO: 3): 

5 1 -GA GGATG CCTTTAT GGATCC A CTCGAG ATCCCAATCCA- 3 ' 
Fokl BamHI Xhol 

20 

This adds the Fok I site which will allow initiation of the sequencing process 
discussed more fully below. 

Tags can be conjugated to cDNAs of existing libraries by standard cloning 
methods. cDNAs are excised from their existing vector, isolated, and then ligated into 

25 a vector containing a repertoire of tags. Preferably, the tag-containing vector is 

linearized by cleaving with two restriction enzymes so that the excised cDNAs can be 
ligated in a predetermined orientation. The concentration of the linearized tag- 
containing vector is in substantial excess over that of the cDNA inserts so that 
ligation provides an inherent sampling of tags. 

30 A general method for exposing the single stranded tag after amplification 

involves digesting a target polynucleotide-containing conjugate with the 5'— >3* 
exonuclease activity of T4 DNA polymerase, or a like enzyme. When used in the 
presence of a single deoxynucleoside triphosphate, such a polymerase will cleave 
nucleotides from 3' recessed ends present on the non-template strand of a double 

35 stranded fragment until a complement of the single deoxynucleoside triphosphate is 
reached on the template strand. When such a nucleotide is reached the 5 f — >3' 
digestion effectively ceases, as the polymerase's extension activity adds nucleotides at 
a higher rate than the excision activity removes nucleotides. Consequently, single 
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12"^: ~ ,ed wi,h n,,cieo,ides - ■»■— * o„,o 

solid phase supports. 

The .echnique may also be used ,o preferentially mesylate imerior Fok I si.es 
of a urge. po ly „ U c,e„,ide while leaving a single Fok I sue a. the .erminus of fe 
polynucleotide unme.hyla.ed. Firs,, me .erminal Fok I si.e is rendered single sttanded 
us ng . polymerase with deoxycy , jdine ^ douHe « d 

1 TnT ,s , n,e,hy ' a,ed ' which ,he sin8le s,randed <~ <• <■ 

«dta a DNA polymerase » the presenee of all four nucleoside triphospha.es, .hereby 
regenerate the Fok I site. Clearly, mis procedure can be generalized ,o 

endonucleases other than Fok I. 

After the oligonucleotide tags are prepared for specific hybridization, e.g. by 
^ndenng them s,„g.e stranded as described above, me P o.ynucleo,ides are mixed 
w,.h mtcropamcles contaming ,he complementary sequences of me tags under 
conditions .ha. favor the formation of perfectly matched duplexes between the tags 
and thetr components. There is extensive guidance in me li.era.ure for creating these 
condmons. Exemplary references providing such guidance include Wetinur Critical 
Renews ,„ Biochemistry and Molecular Biology, 26: 227-259 (1991 ); Sambrook e, 
al Mo,ecu.ar Coning: A Uboratory Manua,, 2nd Edition (Cold Spring Har J 
- Laboratory, New York, ,989,; tmd me like. Preferably, me hybridization conditions 
are suffic.en.ly stnngen. so ma. only perfectly matched sequences form stable 
dup exes. Under such conditions the polynuc.eo.ides specifically hybridized through 
.hen- tags may be l,ga,ed ,o the commentary sequences anached to the 

ZZT l Z Fina " y ' miCropmicl - - ™shed to remove polynucleotides with 
unhgated and/or mismatched tags. 

used ,h W T CP °, miCr0panic,es c <»>ve„,ional.y employed as syndesis supports are 
«d. de„s„y of ug complements on me microparticle surface is W ically greater 
than that necessary for some sequencing operation, That is, in sequencing 
approaches tha, require successive tieaunen, of the attached polynucleotides with a 
vanerv of enzymes, densely spaced polynucleotides may .end to inhibit access of ,he 
reWely bulky enzymes to me polynucleotides. ,„ such cases, the polvnucleo 
are p^ferablv mtxed with ,he micropar.ic.es so ma, tag complement are presen. in 
^mfican, excess, e.g. from ,0:1 ,o 1 00: 1 . or greater, over the polynucleotide^ This 

ZZ> 1 V " ° f P ° ,ynUC,e01idK " - --particle surface „i„ no, be » 
h gh as ,o .nmb,, enzyme access. Preferably, me average in.er-polynuCeo.ide spacmg 
onme m ,cr 0 par.,c,esurfaceis„„ U1 e„rderof 3 0-,OOnm. Guidance in selecting 
ratios forward CPG supports and Ballotini beads ,a type of solid glass support, is 
found ,n Maskos and Southern, Nucleic Acids Research. 20: ,679-1684 ,1992) 
Preferably, for sequencing applications, stiutdard CPC beads of diameter in me range 
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of 20-50 fim are loaded with about 10 5 polynucleotides, and GMA beads of diameter 
in the range of 5-10 |im are loaded with a few tens of thousand of polynucleotides, 
e.g. 4x 10 4 to6x 10 4 . 

In the preferred embodiment, tag complements are synthesized on 

5 microparticles combinatorial^; thus, at the end of the synthesis, one obtains a 

complex mixture of microparticles from which a sample is taken for loading tagged 
polynucleotides. The size of the sample of microparticles will depend on several 
factors, including the size of the repertoire of tag complements, the nature of the 
apparatus for used for observing loaded microparticles--e.g. its capacity, the tolerance 

1 0 for multiple copies of microparticles with the same tag complement (i.e. "bead 
doubles"), and the like. The following table provide guidance regarding 
microparticle sample size, microparticle diameter, and the approximate physical 
dimensions of a packed array of microparticles of various diameters. 

15 

Microparticle diameter 5 10 *im 20 nm 40 jim 

Max. no. 

polynucleotides loaded 

atlperI0 5 sq. 3 x I0 5 1.26 x 10 6 5 x I0 6 

angstrom 

Approx. area of 
monolayer of I0 6 

microparticles .45 x .45 cm I x I cm 2 x 2 cm 4 x 4 cm 

20 The probability that the sample of microparticles contains a given tag complement or 
is present in multiple copies is described by the Poisson distribution, as indicated in 
the following table. 

25 

Table VII 
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Number of 
micropanicJes in 
sample (as fraction 
of repertoire size), 
m 



Fraction of 
repenoire of tag 
complements 
present in 
sample. 



1.000 
.693 
405 
.285 
.223 
.105 
.010 



0.63 
0.50 
0.33 
0.25 
0.20 
0.10 
0.01 



Fraction of 
microparticies in 
sample with unique 
tag complement 
attached, 
m(e- m }/2 

0.37 
0.35 
0.27 
0.21 
0.18 
0.09 
0.01 



Fraction of 
microparticies in 
sample carrying 
same tag 
complement as one 
other microparticle 

in sample 
("bead doubles"). 
m 2 (e" m )/2 

0.18 
0.12 
0.05 
0.03 
0.02 
0.005 



High Specificity Snrting and Panning 

The kinetics of sorting depends on the rate of hybridization of oligonucleotide 
tags to the, tag complements which, in turn, depends on the comp,exity of the tags ,n 
- the hy bnd^anon reaction. Thus, a trade off exists between sorting rate and tag 
complexity, such that an increase in sorting rate may be achieved at the cost of 
reducmg the complexity of the tags involved in the hybridization reaction As 
explained below, the effects of this trade off may be ameliorated by "panning - 

Specificity of the hybridizations may be increased by taking a sufficiently 
small sample so that both a high percentage of tags in the sample are unique and the 
nearest ne.ghbors ofsubstantia.lv all the tags in a sample differ bv at least two words 
Th.s latter condition may be met by taking a sample that contains a number of tag- ' 
po.ynucleot.de conjugates that is about 0. . percent or less of the size of the repenoire 
bong emp.oyed. For example, if tags are constructed with eight words selected from 
Table II, a repertoire of 8 8 . or about . .67 x . 0 7 tags and tag complements are 
produced. In a library of tag-cDNA conjugates as described above, a 0 1 percent 
sample means that about ,6,700 different tags are present. If , hlS were loaded duec.lv 
onto a repertoire-equivalent of microparticies. or in this example a sample of 1 67 x ' 
0 m.cropart.cles, then only a sparse subset of the sampled microparticies would be 
loaded. The densuy of .oaded micropart.cles can be .ncrease-for example, for more 
efficient sequencing-by undertaking a "panning" step in which the sampled tag- 
cDNA conjugates are used to separate loaded microparticies from unloaded 
rmcroparticles. Thus, in the example above, even though a "0. 1 percent" sample 
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contains only 1 6,700 cDNAs, the sampling and panning steps may be repeated until 
as many loaded microparticles as desired are accumulated. 

A panning step may be implemented by providing a sample of tag-cDNA 
conjugates each of which contains a capture moiety at an end opposite, or distal to, 
5 the oligonucleotide tag. Preferably, the capture moiety is of a type which can be 
released from the tag-cDNA conjugates, so that the tag-cDNA conjugates can be 
sequenced with a single-base sequencing method. Such moieties may comprise 
biotin, digoxigenin, or like ligands, a triplex binding region, or the like. Preferably, 
such a capture moiety comprises a biotin component. Biotin may be attached to tag- 

1 0 cDNA conjugates by a number of standard techniques. If appropriate adapters 

containing PCR primer binding sites are attached to tag-cDNA conjugates, biotin may 
be attached by using a biotinylated primer in an amplification after sampling. 
Alternatively, if the tag-cDNA conjugates are inserts of cloning vectors, biotin may be 
attached after excising the tag-cDNA conjugates by digestion with an appropriate 

1 5 restriction enzyme followed by isolation and filling in a protruding strand distal to the 
tags with a DNA polymerase in the presence of biotinylated uridine triphosphate. 

After a tag-cDNA conjugate is captured, it may be released from the biotin 
moiety in a number of ways, such as by a chemical linkage that is cleaved by 
-reduction, e.g. Herman et al, Anal. Biochem., 156: 48-55 (1986), or that is cleaved 

20 photochemical ly, e.g. Olejnik et al, Nucleic Acids Research, 24: 361-366 (1996), or 
that is cleaved enzymatically by introducing a restriction site in the PCR primer. The 
latter embodiment can be exemplified by considering the library of tag-polynucleotide 
conjugates described above: 

25 S'-RCGACCAICjW^WlgGGmxg- cDNA -NNNR 

GGT [G, W, W, W] 9 CC[A] 19 - r DNA -NNNYCTAG-5 1 

The following adapters may be ligated to the ends of these fragments to permit 
amplification by PCR: 

30 

5 ■ - xxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXXXXYGAT 
35 Right Adapter 



GATCZZACTAGTZZZZZZZZZZZZ-3 r 
40 ZZTGATCAZZZZZZZZZZZZ 
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Left Adapter 
ZZTGATCAZZZZZZZZZZZZ-5 • -biotin 
Left Primer 

where -ACTAGT" is a Spe . recognition site (which lMves . 

Z L , Se " UH,Cin8) - "* X ' S ^ Z ' S « - to, 

the an„eal,„g and dtssociation temperatures of the respective primers are 

approxunarely ,he same. After ligation of the adapters and amplification by PCR 

usmg the btotinylated primer, the tags of ft. conjugates are rendered singie stianded 

by the exonuclease activity of T4 DNA polymerase and conjugates are combined with 

TTh* "I r Cr0Pa " ideS - C ' 8 3 repenoire ">""^- w»» <aE comp.emen,s 
attached. After annealing under sttingen, conditions (,o m in imize m is. atl achment of 

the conjugates are preferably ligated to their tag com pl eme„.s and the loaded 
mtcropamcles are separated from the unloaded microparticles by capture with 
avid.nated magnetic beads, or like capture technique. 

, „ <nn f^TTf IO ** eXaroPle ' °" S Pr0CeSS reSU " S in ,he accumulation of about 

I .500 (-.6,700 x .63, loaded microparticles with different tags, which may be 

o e ^ir t maene,ic ^ by ciravage wi,h * 

II T *"* *" ° f "*"*"■«'*• — ^-cDNA conjugates. 4-5 x 105 
cDNAs can be accumulated by pooling the released microparticles. The pooled 

25 tTh^uT'" ^ ' hCn SimU " aneOUS ^ "*»««" a single-base sequencing 

Determining how many Umes to repeat the sampling and panning s,e P s-or 
more generally, determining how many cDNAs to analyze, depends on one's 
objective. If the objective is to monitor the changes to abundance of relatively 
common sequences, e.g. making „ P 5% or more of a population, then relatively small 
samples. ,.e. a smal, fraction of the total population size, may allow statistical, 
s.gn.flcan, estimates of relative abundances. On the other hand, if one seeks to 
momtor the abundances of rare sequences, eg making up 0..% or less of a 
population, then large samples are required. Generally, there is a direct relationship 

otrz;: r ^ ,he reiiabiii,y of ^ ~ ° { '«*™ — — 

on the sample. There ,s extenstve guidance in the literature on determining 
appropriate sample sizes for making reliable statistical estimates, e.g. Roller et al 
Nucletc Acids Research, 23:185-191 (1994): Good, Biometrika. 40: 16-264 (1953, 
Bunge e, al. J. Am. Stat. Assoc., 88: 364-373 (1993); and the like. Preferably fo, 
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monitoring changes in gene expression based on the analysis of a series of cDNA 
libraries containing 10 5 to 10 8 independent clones of 3.0-3.5 x 10 4 different 
sequences, a sample of at least 10 4 sequences are accumulated for analysis of each 
library. More preferably, a sample of at least 10 5 sequences are accumulated for the 
5 analysis of each library; and most preferably, a sample of at least 5 x 10 5 sequences 
are accumulated for the analysis of each library. Alternatively, the number of 
sequences sampled is preferably sufficient to estimate the relative abundance of a 
sequence present at a frequency within the range of 0.1% to 5% with a 95% 
confidence limit no larger than 0.1% of the population size. 

10 

Single Base DNA Sequencing 
The present invention can be employed with conventional methods of DNA 
sequencing, e.g. as disclosed by Hultman et ah Nucleic Acids Research, 17: 4937- 
4946 (1989). However, for parallel, or simultaneous, sequencing of multiple 

1 5 polynucleotides, a DNA sequencing methodology is preferred that requires neither 
electrophoretic separation of closely sized DNA fragments nor analysis of cleaved 
nucleotides by a separate analytical procedure, as in peptide sequencing. Preferably, 
the methodology permits the stepwise identification of nucleotides, usually one at a 
- time, in a sequence through successive cycles of treatment and detection. Such 

20 methodologies are referred to herein as "single base" sequencing methods. Single 
base approaches are disclosed in the following references: Cheeseman, U.S. patent 
5,302,509; Tsien et al, International application WO 91/06678; Rosenthal et al, 
International application WO 93/21340; Canard et al, Gene, 148: 1-6 (1994); and 
Metzker et al, Nucleic Acids Research, 22: 4259-4267 (1994). 

25 A "single base" method of DNA sequencing which is suitable for use with the 

present invention and which requires no electrophoretic separation of DNA fragments 
is described in International application PCT/US95/03678. Briefly, the method 
comprises the following steps: (a) ligating a probe to an end of the polynucleotide 
having a protruding strand to form a ligated complex, the probe having a 

30 complementary protruding strand to that of the polynucleotide and the probe having a 
nuclease recognition site; (b) removing unligated probe from the ligated complex; (c) 
identifying one or more nucleotides in the protruding strand of the polynucleotide by 
the identity of the ligated probe; (d) cleaving the ligated complex with a nuclease; and 
(e) repeating steps (a) through (d) until the nucleotide sequence of the polynucleotide. 

35 or a portion thereof, is determined. 

A single signal generating moiety, such as a single fluorescent dye, may be 
employed when sequencing several different target polynucleotides attached to 
different spatially addressable solid phase supports, such as fixed microparticles, in a 
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parallel sequencing opera,io„. This may be accomplished by providjng fou , ^ 
probes that are applied sequential,, lo of ^rge, polynucj,,^ 

d,ffere„, nncroparUc.es. An exempt se, of such probes are shown belo" 



Set 1 



Set 2 



Set 3 



Set 4 



ANNNN . . . NN 

N. . . NNTT . . . T* 

dCNNNN . . . NN 

N. . .NNTT. . . T 

dGNNNN . . . NN 

N . . . NNTT . . . T 

dTNNNN . . . NN 

N . . . NNTT . . . T 



NNTT . . . T »...MNTT...T N . . .NNTT. . . T 

CNNNN . . . NN wno „ ilM 

N...NNTT T* dCN ^--- NN dCNNNN . . . NN 

NNTT . . . T N...NNTT...T N...NMTT...T 

dGNNNN . . . NN GNNNN . . . NN rirutJUM 

N . . NNTT T m dGNNNN... NN 

••■KW11...T N. . . NNTT T« 

N...NNTT...T 

dTNNNN . . . NN htuhmm 

N NNTT t * • • TNNNN... NN 

N . . . NNTT . . . T N...NNTT...T N . . . NNTT . . . T * 



positions in the protruding strand are fiii^ k , 0ther 

Preferably, such 3'-terminal nucleotides are dideoxynucleotides In J, u t 

identified bv the signa. generated bvTe 1 h I T """"S " A " m 

cleavage. the probeTof " aD l> d l^" " ^ ™^ - 

H 8 a.ed complL wi, h ^^^^.7%''?^^ f<mnin8 
the probes of sets 3 „nH a ™ ma ""6 '■> C are tdenltfied by location. Similarly. 

probes of sets 3 and 4 are apphed and locations of positiye signals identify ThiV 
process o sequentially applying the four sets of probes continues J, t e^t. 
number of nucleot.des are identified on the urge, polynucleotides. Clejyone^f 
or tnarysk,,, could construe similar sets of probes ,ha, co„,d have m^vXl 
uc as navtng protrudtng strands of different lengths, different moieties to l 
hgatton of unlabeled probes, different means for labeling probes, and the like 
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A pparatus for Sequencing Populations of Polynucleotides 
An objective of the invention is to sort identical molecules, particularly 
polynucleotides, onto the surfaces of microparticles by the specific hybridization of 
tags and their complements. Once such sorting has taken place, the presence of the 

5 molecules or operations performed on them can be detected in a number of ways 
depending on the nature of the tagged molecule, whether microparticles are detected 
separately or in "batches," whether repeated measurements are desired, and the like. 
Typically, the sorted molecules are exposed to ligands for binding, e.g. in drug 
development, or are subjected chemical or enzymatic processes, e.g. in polynucleotide 

1 0 sequencing. In both of these uses it is often desirable to simultaneously observe 

signals corresponding to such events or processes on large numbers of microparticles. 
Microparticles carrying sorted molecules (referred to herein as "loaded" 
microparticles) lend themselves to such large scale parallel operations, e.g. as 
demonstrated by Lam et al (cited above). 

1 5 Preferably, whenever light-generating signals, e.g. chemiluminescent, 

fluorescent, or the like, are employed to detect events or processes, loaded 
microparticles are spread on a planar substrate, e.g. a glass slide, for examination with 
a scanning system, such as described in International patent applications 
• PCT/US91/09217,PCT/NL90/00081,andPCT/US95/01886. The scanning system 

20 should be able to reproducibly scan the substrate and to define the positions of each 
microparticle in a predetermined region by way of a coordinate system. In 
polynucleotide sequencing applications, it is important that the positional 
identification of microparticles be repeatable in successive scan steps. 

Such scanning systems may be constructed from commercially available 

25 components, e.g. x-y translation table controlled by a digital computer used with a 
detection system comprising one or more photomultiplier tubes, or alternatively, a 
CCD array, and appropriate optics, e.g. for exciting, collecting, and sorting 
fluorescent signals. In some embodiments a confocal optical system may be 
desirable. An exemplary scanning system suitable for use in four-color sequencing is 

30 illustrated diagrammatically in Figure 5. Substrate 300, e.g. a microscope slide with 
fixed microparticles, is placed on x-y translation table 302, which is connected to and 
controlled by an appropriately programmed digital computer 304 which may be any of 
a variety of commercially available personal computers, e.g. 486-based machines or 
PowerPC model 7100 or 8100 available form Apple Computer (Cupertino, CA). 

35 Computer software for table translation and data collection functions can be provided 
by commercially available laboratory software, such as Lab Windows, available from 
National Instruments. 
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Substrate 300 and table 302 are operational* associated with microscope 306 
hav mg one „ objectjve lense$ whjch ^ 306 

3e 3 - 300 E «iu. ion bea^ 0 fr0m 

l.ght source 312, wh,ch ,s preferably a laser, is directed to beam splitter 314 e „ . 

£T "f ^ 1,16 * hr ° Ueh ™>* ^ - ^ e ens 

filters, fiber opt,cs, rotatmg mirrors, fixed position mirrors and lenses diffraction 

gratings, and the like. As illustrated in Fieure 1 a- . u dlffractl °" 

fluorescence 316 tn W L ? ^ distribution optics 3 1 8 directs 

orescence 3 16 to four separate photomultiplier tubes. 330, 332 334 and 336 
who S 1S tn directed to ^ ^ photon countJ * ^ and 3* 

356. The output of the photon counters is collected by computer 304 wh er , u 
stored, analyzed, and viewed on video 360 Alters 1 7 \ * 
3 1 8 could be a Hiffr *• Alternately, signal distribution optics 

^ ^ 8 could be a diffraction grating which directs fluorescent signal 3 1 8 onto a CCD 

d^S^t^ of ; he positionai ,oca,ization in ~ ** Wi » 

, aiarg e extent, the resolution for separating closely spaced 
rn^pamce, Prc ferab, y> the scanning sy s.em S capab Xso,™ 

between a out ten thousand to one hundred thousand micropLcles ^ m 2 

In sequencing applications, loaded microparticles can be fix J L * ' * 
of a substrate in variety of wav, Th.fi u ° the SUrface 

micropartides to M Je^ E o su^,^ " " " ^ '° a " OW "* 
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in turn may be cross-linked to avidin, again using conventional chemistries, to form 
an avidinaied surface. Biotin moieties can be introduced to the loaded microparticles 
in a number of ways. For example, a fraction, e.g. 10-15 percent, of the cloning 
vectors used to attach tags to polynucleotides are engineered to contain a unique 
5 restriction site (providing sticky ends on digestion) immediately adjacent to the 

polynucleotide insert at an end of the polynucleotide opposite of the tag. The site is 
excised with the polynucleotide and tag for loading onto microparticles. After 
loading, about 10-15 percent of the loaded polynucleotides will possess the unique 
restriction site distal from the microparticle surface. After digestion with the 

1 0 associated restriction endonuclease, an appropriate double stranded adaptor 

containing a biotin moiety is ligated to the sticky end. The resulting microparticles 
are then spread on the avidinated glass surface where they become fixed via the 
biotin-avidin linkages. 

Alternatively and preferably when sequencing by ligation is employed, in the 

1 5 initial ligation step a mixture of probes is applied to the loaded microparticle: a 

fraction of the probes contain a type lis restriction recognition site, as required by the 
sequencing method, and a fraction of the probes have no such recognition site, but 
instead contain a biotin moiety at its non-ligating end. Preferably, the mixture 
* comprises about 10-15 percent of the biotinylated probe. 

20 In still another alternative, when DNA-loaded microparticles are applied to a 

glass substrate, the DNA may nonspecifically adsorb to the glass surface upon several 
hours, e.g. 24 hours, incubation to create a bond sufficiently strong to permit repeated 
exposures to reagents and washes without significant loss of microparticles. 
Preferably, such a glass substrate is a flow cell, which may comprise a channel etched 

25 in a glass slide. Preferably, such a channel is closed so that fluids may be pumped 
through it and has a depth sufficiently close to the diameter of the microparticles so 
that a monolayer of microparticles is trapped within a defined observation region. 

Identification of Novel Polynucleotides 
30 in cDNA Libraries 

Novel polynucleotides in a cDNA library can be identified by constructing a 
library of cDNA molecules attached to microparticles, as described above. A large 
fraction of the library, or even the entire library, can then be partially sequenced in 
parallel. After isolation of mRNA, and perhaps normalization of the population as 
35 taught by Soares et al, Proc. Natl. Acad. Sci. ? 91 : 9228-9232 (1994), or like 

references, the following primer may by hybridized to the polyA tails for first strand 
synthesis with a reverse transcriptase using conventional protocols (SEQ ID NO: 1 ): 
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5 ' -mRNA- fA] n -3' 

rT] is- [primer site) -GG [W, W, W, C] 9 ACCAGCTGATC-5 • 

where [W, W,W,C]o represents a tag as described above " ACCAGCTC A TP" • 
> "-nee ^^^^^^^^ 
a sequence common to all members of the library that is ,ater used as a priZ 
bmdmg sate for amplifying polynucleotides of interest by PCR 

After reverse transcription and second strand synthesis by conventional 
echoes the double stranded fragments are inserted into a cloning vectr " 

3^"^^*^ Preferably, as desridTove the 
clonrng ector , s constructed so that the tagged cDNAs can be excised with an ' 
endonuclease, such as Fok I that will allow i mm ~r ♦ 

single base methoH * immediate sequencmg by the preferred 

single base method after sortmg and ligation to microparticles 

Stepwise sequencing is then carried out simultaneously on the whole lib™™ 
or one or more large fractions of the library, in accordance wi^he ^' 
suffice number of nucleotides are .de^'ed 0.^^^^ ^ ' 

»fIT; he "** ° r8aniSm fr0m ^ * derived. For 

example, f the hbrary » denved from mammalian mRNA then a randomly selected 
sequence 14-15 nucleotides long is exoected to h» w • ^""V selected 

derive <■ i suniclcm,or ™q"c representation in a library 

der ved from ba«e„a. or Cher !„wer organisms. PreferaWy. a. least 20-30 

compared ,o known sequences to identify uniaue cDNAs 

Unrque cDNAs are then isolated by conventional technioues e 8 construct^ 
a probe from the PCR amplicon produced with primers directed ,o * '„ ""^"J 
*e portion of the cDNA whose sequence was l^J^Z^ZT 
— .0 ,de„,«v the cDNA in a Horary using a conventional screening pro „t„ " 

The above method for identifying new cDNAs rnav ak n h TT 
mRNA populations, enher m lsol a,ed ^»^.7^^r 8en>rin ' 
dynam.cally changing P o P „ la ,io„. Partial sequence information is ob "led 
sunuhaneously from a rarge sanpre, e.g. ten , 0 a hundred thousand o^ of 
CDNAs anached ,0 separate microparticles as described in ,he above mZd 
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Example 1 

Construction of a Tae Library 
An exemplary tag library is constructed as follows to form the chemically 
5 synthesized 9-word tags of nucleotides A, G, and T defined by the formula: 

3'-TGGC-[ 4 <A,G,T) 9 ]-CCCCp 

where "[^AjGJ)^" indicates a tag mixture where each tag consists of nine 4-mer 
10 words of A, G, and T; and M p" indicate a 5' phosphate. This mixture is ligated to the 
following right and left primer binding regions (SEQ ID NO: 4 and SEQ ID NO 5): 



15 



5 • - AGTGGCTGGGCATCGGACCG 5 1 - GGGGCCCAGTCAGCGTCGAT 

TCACCGACCCGTAGCCp GGGTCAGTCGCAGCTA 

LEFT RIGHT 



The right and left primer binding regions are ligated to the above tag mixture, after 
which the single stranded portion of the ligated structure is filled with DNA 
20 'polymerase then mixed with the right and left primers indicated below and amplified 
to give a tag library (SEQ ID NO: 6). 



Left Pruner 

25 

5 1 - AGTGGCTGGGCATCGGACCG 

5 ' - AGTGGCTGGGCATCGGACCG- [ 4 (A, G, T) 9] -GGGGCCCAGTCAGCGTCGAT 
TCACCGACCCGTA GCCTGGC - [ 4 (A, G, T) 9 ] -CCCCGGGTCAGTCGCAGCTA 

30 

CCCCGGGTCAGTCGCAGCTA- 5 ' 
Right Primer 



35 The underlined portion of the left primer binding region indicates a Rsr II recognition 
site. The left-most underlined region of the right primer binding region indicates 
recognition sites for Bsp 1201, Apa L and Eco O 1091, and a cleavage site for Hga I. 
The right-most underlined region of the right primer binding region indicates the 
recognition site for Hga I. Optionally, the right or left primers may be synthesized 

40 with a biotin attached (using conventional reagents, e.g. available from Clontech 
Laboratories, Palo Alto, CA) to facilitate purification after amplification and/or 
cleavage. 
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primer binding site Ppu Ml site 




-CAAATTTG- CCTAGG- AGAAGGAGAAGGAGAAGG- 
t t 

Bam HI site 

Pme I site 



15 

The plasmid is cleaved with Ppu MI and Pme 1 (to give a Rsr H-compatible end and a 
flush end so that the insert is oriented) and then methylated with DAM methylase. 
The tag-containing construct is cleaved with Rsr II and then ligated to the open 
plasmid, after which the conjugate is cleaved with Mbo I and Bam HI to permit 
20 ligation and closing of the plasmid. The plasmid is then amplified and isolated and 
* used in accordance with the invention. 



Example 3 

Changes in Gene Expression Profiles in Liver Tissue of Rats 

25 Exposed to Various Xenobiotic Agents 

In this experiment, to test the capability of the method of the invention to 
detect genes induced as a result of exposure to xenobiotic compounds, the gene 
expression profile of rat liver tissue is examined following administration of several 
compounds known to induce the expression of cytochrome P-450 isoenzymes. The 

30 results obtained from the method of the invention are compared to results obtained 

from reverse transcriptase PCR measurements and immunochemical measurements of 
the cytochrome P-450 isoenzymes. Protocols and materials for the latter assays are 
described in Morris et al, Biochemical Pharmacology , 52: 781-792 (1996). 

Male Sprague-Dawley rats between the ages of 6 and 8 weeks and weighing 

35 200-300 g are used, and food and water are available to the animals ad lib. Test 
compounds are phenobarbital (PB), metyrapone (MET), dexamethasone (DEX), 
clofibrate (CLO), corn oil (CO), and P-naphthoflavone (BNF), and are available from 
Sigma Chemical Co. (St. Louis, MO). Antibodies against specific P-450 enzymes are 
available from the following sources: rabbit anti-rat CYP3A1 from Human Biologies, 

40 Inc. (Phoenix, AZ); goat anti-rat CYP4A1 from Daiichi Pure Chemicals Co. (Tokyo, 
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Japan); monoclonal mouse ami-rat CYP1 AI. monoclonal mouse anti-ntt CYP2C1 1 
goa,an,,ra,CYP2El, and monoclonal mouse an,i- ra , CYP2B1 from Oxford 
Biochem ical Research, !nc. (Oxford, MI). Secondary antibodies (goa, anti-rabbi, Ie0 
rabb,, ami-goa, IgG and goa, ami-mous, IgG) are available from Jackson ' 
= ImmunoResearch Uboratories (West Grove, PA). 

,,„„ ^ imalSaread ' ninis,ercd ' :iU,erPB ( l0( >"'g^g),BNF(10()mg*g) MET 
(.00 mg*g>, DEX (100 mgfcg), or CLO (250 mg/kg, fo, 4 consecutive^ 
m^pentonea, injection following . dosing regimen similar ,o ma, desert 

H 2 0 and CO are used as controls. Two hours following me las, injec.ion (day 4) 

memod? 131 ^ USing a —«*■«« of me 

5 Too rr'f^ * " Bi °' eChni,UeS ' ' ' : 32 " 27 Approximately 
£200 mg o hver tissue ,s homogenized in me RNA extraction buffer descri bed bv 
Xie e, a, ,„ , solaK tota , ^ ^ ^ ^ ^ > 

d^ylpy r ocarb„ n a K . tt ea,ed wa,er, auantif.ed spec,ropl,o,ome«rically a, 260 „m and 
adjus,ed ,o a co„cenua,ion ofl 00 ug/ml. Toul RNA is slored in 

- dielhylpyrocarbonale-lreated water for un to I v„«r „. mor- -.t. 
J degradation rtpt, j ror UD 10 1 at -70o C without any apparent 

degradaMon RT-PCR and sequencing are performed on samples from these 

preparations. 

For sequencing, samples of RNA corresponding ,o about 0.5 „g of poly,A)+ 
RN A are used to construe, libraries of tag-cDNA conjugates following the proCi 
d scribed ,n the section e„,i„ed "Anaching Tags ,o Polynucleotides for Sorting™ 
Solid Phase Supports,- with me following exception: the tag repertoire is consJTed 
from six 4-nucleot de words from Table II Th„. .t . «>ns,ruc,w 
8' or about 2 6, ins p u complexity of me repertoire is 

8 orabou,2.6x .05. For each tag-cDN A conjugate library constructed, ten samples 
of afciu, ttn thousand clones are uken for amplincanon and sorting. Each of uT 
amplified samples is separately applied to a fixed monolayer of abou, .0* ,0 1 
d,ame,er GMA beads conuining u,g complements. Tha, is, me -sample" of uT 

Z T " P ° PUlati0n °" "* * four Z me 

ota, ,ze of the repertoire, Uius ensuring mere ,s a high probabi,i,y tha, each o t T 

sampled tag-cDNA conjugates will f,„d its tag complement on the monolayer 

he oligonucleotide tags of the amplified samples are rendered single sanded as 

described above, the lag-cDNA onjuga.es of the samples are separately appll the 

monolayers under conditions ma, permit specific hybridization only between 

oligonucleotide tags and tag complements forming perfectly matched duplexes 

Concentrates of me amplified samples and hybridization times are selected to 
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permit the loading of about 5 x 10 4 to 2 x 10 5 tag-cDNA conjugates on each bead 
where perfect matches occur. After ligation, 9-12 nucleotide portions of the attached 
cDN As are determined in parallel by the single base sequencing technique described 
by Brenner in International patent application PCT/US95/03678. Frequency 
5 distributions for the gene expression profiles are assembled from the sequence 
information obtained from each of the ten samples. 

RT-PCRs of selected mRNAs corresponding to cytochrome P-450 genes and 
the constitutively expressed cyclophilin gene are carried out as described in Morris et 
al (cited above). Briefly, a 20 fiL reaction mixture is prepared containing lx reverse 

10 transcriptase buffer (Gibco BRL), 10 nM dithiothreitol, 0.5 nM dNTPs, 2.5 jaM oligo 
d(T)| 5 primer, 40 units RNasin (Promega, Madison, Wl), 200 units RNase H-reverse 
transcriptase (Gibco BRL) ¥ and 400 ng of total RNA (in diethylpyrocarbonate-treated 
water). The reaction is incubated for 1 hour at 37°C followed by inactivation of the 
enzyme at 95°C for 5 min. The resulting cDNA is stored at -20°C until used. For 

15 PCR amplification of cDNA, a 10 |iL reaction mixture is prepared containing lOx 
polymerase reaction buffer, 2 mM MgCl 2 > 1 unit Taq DNA polymerase (Perkin- 
Elmer, NorwaIk ? CT), 20 ng cDNA, and 200 nM concentration of the 5' and 3* 
specific PCR primers of the sequences described in Morris et al (cited above). PCRs 
-are carried out in a Perkin-Elmer 9600 thermal cycler for 23 cycles using melting, 

20 annealing, and extension conditions of 94°C for 30 sec, 56°C for 1 min., and 72°C 
for 1 min., respectively. Amplified cDNA products are separated by PAGE using 5% 
native gels. Bands are detected by staining with ethidium bromide. 

Western blots of the liver proteins are carried out using standard protocols 
after separation by SDS-PAGE. Briefly, proteins are separated on 10% SDS-PAGE 

25 gels under reducing conditions and immunoblotted for detection of P-450 isoenzymes 
using a modification of the methods described in Harris et al, Proc. Natl. Acad. ScL 
88: 1407-1410 (1991). Protein are loaded at 50 ng/lane and resolved under constant 
current (250 V) for approximately 4 hours at 2°C. Proteins are transferred to 
nitrocellulose membranes (Bio-Rad, Hercules, CA) in 1 5 mM Tris buffer containing 

30 120 mM glycine and 20% (v/v) methanol. The nitrocellulose membranes are blocked 
with 2.5% BSA and immunoblotted for P-450 isoenzymes using primary monoclonal 
and polyclonal antibodies and secondary alkaline phosphatase conjugated anti-IgG. 
Immunoblots are developed with the Bio-Rad alkaline phosphatase substrate kit. 
The three types of measurements of P-450 isoenzyme induction showed 

35 substantial agreement. 



-43- 



BNSDOCID: cWO 9713B77A1 I > 



WO 97/13877 



PCTYUS96/16342 



Program minxh 
c 



APPENDIX la 

Exemplary ^m p ,^ program for f^n^,,;^ 
minimally r r OS s hvhrirfiyinp r cm c 
(single stranded tag/single stranded tag complement) 



C 

c 



c 
c 



c 
c 



c 
c 



c 



900 
c 



d?^n! r * 2 S " bl(6) ' mSetl ( 10 00,6), m set2( 10 00,6, 

dimension nbase(6) ^,o> 



writer,*) 'ENTER SUBUNIT LENGTH ' 
read(*, 100)nsub 
100 format (il ) 



open ( 1 , f il e = • sub4 . dat . , f orn)= . f orn|at ted , ^ , 



nset=0 

do 7000 ml = l, 3 

do 7000 m2=l, 3 

do 7000 m3=l, 3 

do 7000 m4=l, 3 
subl (l)»ml 
subl (2)=m2 
subl [ 3)=m3 
subl (4)=ro4 



ndiff =3 



Generate set of subunits differing frc 
c subl by at least ndiff nucleotides 

Save in msetl . 



jj-1 

do 900 j«l,nsub 

msetl (1, j ) =subl { j ) 



do 1000 kl-1, 3 
do 1000 k2=l, 3 

do 1000 k3 = l, 3 

do 1000 k4 = l,3 



nbase (1 ) =kl 
nbase (2)=k2 
nbase(3)=k3 
nbase ( 4 ) =k4 



-44- 



WO 97/13877 



PCT/US96/16342 



1200 

c 
c 



n=0 

do 1200 j=l,nsub 
if (subl { j ) .eq 
1 subl ( j ) . eq 

3 subf ( j } .eq 

n=n+ 1 
endif 
continue 



if (n . ge . ndif f ) then 



. and . 
. and . 
. and . 



nbase ( j ) 
nbase ( j ) 
nbase ( j ) 



ne 
ne 
ne 



1 
2 
3) 



or . 
or . 
then 



c 
c 

c 
c 
c 
c 



1100 

c 
c 

1000 

c 

c 

1325 



If number of mismatches 
is greater than or equal 
to ndiff then record 
subunit in matrix mset 



j j-j j + 1 
do 1100 i=l,nsub 

mset 1 ( j j , i ) =nbase ( i ) 
endif 



cont inue 



do 1325 j2=l,nsub 

mset 2 (1, j2) =msetl { 1, j2) 

mset 2 (2, j2> =mset 1 (2, j2) 



c 
c 

c 
c 
c 
c 
c 
c 
c 
c 



Compare subunit 2 from 
msetl with each successive 
subunit in msetl, i.e. 3, 
4,5, ... ere. Save those 
with mismatches .ge. ndiff 
in matrix mset2 starting at 
position 2. 

Next transfer contents 
of mset2 into msetl and 
start 

comparisons again this time 
starting with subunit 3. 
Continue until all subunits 
undergo the comparisons. 



c 
c 

1700 



npass=0 



continue 
kk-npass+2 
npass=npass+ 1 
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do 



2 
2 



1600 



1625 



1500 



1500 m*=npass + 2, j j 

do 1600 j=l,nsub 

if(msetl(npass + l /j ). eq .i. and . Insetlfm 

; s «j;npass + l f3 ,.eq.a.and.msetl (m# j 
msetl (npass+1, j , . eq . 3 . and . mset 1 (m, j ) 

endif 
continue 
if (n.ge.ndiff ) then 
kk=kk+l 

do 1625 i~l,nsub 

mset2 fkk, i) -msetl (m r i) 

endif 
continue 



ne 
ne 



or . 
or . 



ne.3) then 



c 
c 
c 
c 
c 
c 



2000 



7009 

7008 
7010 



120 
7000 

c 

c 



*k is the number of subunits 
stored in mset2 

Transfer contents of mset2 
into msetl for next pass. 



do 2000 k=l,kk 

do 2000 m=l,nsub 

msetl f k,m)=mset2 ( k,m) 
if (kk.lt. jj) then 
jj=kk 
goto 1700 
endif 



nset =nset+] 
write {1, 7009) 

format { / ) 
do 7008 k=l, kk 

write (*, *) 

writer, 120) kk,nset 

'°™"LY SUbU "' tS N„-. i5l 

close { 1 ) 



end 



r * * * * W * 1 



r * + * * * + + , 
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APPENDIX lb 

Exemplary computer program for generating 
minimally cross hybridizing sets 
(single stranded tag/single stranded tag complement) 



Program tagN 

c 

c 

c Program tagN generates minimally cross-hybridizing 

c sets of subunits given i) N--subunit length, and ii) 

c an initial subunit sequence. tagN assumes that only 

c 3 of the four natural nucleotides are used in the tags. 

c 
c 

character*! subl (20) 

integer*2 mset ( 10000 , 20 ) , nbase(20) 

c 
c 

write (*,*) * ENTER SUBUNIT LENGTH 1 

read ( * , 100) nsub 
100 format (i2) 

c 
c 

write (*,*) 'ENTER SUBUNIT SEQUENCE 1 

read(*, 110) { subl ( k) , k=l , nsub) 
110 format (20al) 

c 
c 

ndif f =10 

c 

c Let a=l c=2 g=3 & t=4 



800 



do 800 kk=l , nsub 

if (subl (kk) .eq. 'a • ) then 

mset ( 1 , kk) =1 

endi f 

if (subl (kk) .eq. ' c' ) then 
mset (1 , kk) =2 
endi f 

i f ( subl ( kk } . eq . ' g 1 ) then 
mset (1, kk) =3 
endi f 

if (subl { kk) .eq. 1 t • ) then 
mset ( 1 , kk ) ~4 
endi f 

continue 



Generate set of subunits differing from 
subl by at least ndiff nucleotides. 



j j = l 

c 

do 1000 k i = 1 , 3 



-47- 



Q7ir?P77A1 i 



WO 97/13877 

PCT/US96/16342 

do 1000 k2=l, 3 
do 1000 k3=l, 3 
do 1000 k4=l , 3 
do 1000 k5=l, 3 
do 1000 k6=l, 3 
do 1000 k7 = i, 3 
do 1000 k8=l,3 
do 1000 k9=l,3 
do 1000 kl0=l ? 
do 1000 kll«l,3 ' 
do 1000 kl2-l,3 
do 1000 kl3=l, 3 
do 1000 kl4=l,3 
do 1000 kl5=l,3 
do 1000 kl6=l, 3 
dc 1000 Jcl7 = l,3 
do 1000 kl8=l,3 
do 1000 kl9=l, 3 
do 1000 k20-l,3 



c 
c 



nbase <l)=kl 
nbase (2) =k2 
nbase f 3) =k3 
nbase { 4 ) =k4 
nbase (5) =k5 
nbase (6) =k6 
nbase (7j =k7 
nbase (8) =k8 
nbase (9) = k9 
nbase (10 j =kl0 
nbase (11) =kll 
nbase ( 12) =kl2 
nbase(13)-kl3 
nbase (14) =kl4 
nbase (15) -kl5 
nbase(16)=kl6 
nbase (17)=kl7 
nbase (18)=kl8 
nbase (19) = k!9 
nbase (20) = k20 



do 1250 nn = l, j j 



1200 



n=0 

do 1200 j=l,nsub 

if (mset (nn, j ) .eq 
mset (nn, j ) .eq 
mset (nn, j ) 
mset (nn, j ) 
n=n + 1 
endif 
continue 



1 
2 

eq. 3 
eq. 4 



-and. nbase ( j ) . ne. 1 
. and. nbase ( j ) . ne. 2 
• and. nbase ( j ) . ne . 3 
and. nbase(j) .ne.4) 



. or . 
. or . 
. or . 
then 



12i 

c 



iftn.lt. 
goto 
endif 

continue 



ndif f ) 
1000 



then 



j+1 

writer, 130) (nbase(i) ,i = 
do 1100 i»l,nsub 



1, nsub) , i j 
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mset ( j j , i ) =nbase ( i ) 
HOC continue 



100C continue 
c 

write (*, * } 
130 format (10x,20(lx,il),5x,i5) 

write ( * , * ) 

write(*,12C) jj 
120 format (lx, 'Number of words= , ,i5J 

c 



c 
c 
c 



end 
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c 
c 



APPENDIX Ic 

Exemplary computer p ro ^m f or p enPrat ,„ c 
minimally rmcc hvbririmn p c»tc 
(double stranded tag/single stranded tag complement) 

Program 3tagN 



characters subl (20) 

integers msec (10000, 20) , nbase(20) 



write (*, *) 'ENTER SUBUNIT LENGTH ' 
read <*.100)nsub 
100 format (i2) 



write ( * , * ) ■ ENTER SUBUNIT SEQUENCE a * „ n „, . 
readC. 110) (subl <k> , k-1. nsub) 9 ° nly 

format f?n a i \ ; 



*10 format (20al) 

c 

ndiff=lo 

c 

^ Let a=l and g=2 

do 800 kk=l,nsub 

if (subl(kk) , eq . 'a') then 

mset (1, kk)-l 

endif 

if (subl(kk) .eq. -g.) then 

™set (1 , kk) =2 

ft no * ■ endif 

frt'O continue 



33-1 

do 1000 kl=l,3 
do 1000 k2=l, 3 
do 1000 k3=l, 3 
do 1000 k4-l, 3 
do 1000 kS=i, 3 
do 1000 k6-l, 3 
do 1000 k7=a, 3 
do 1000 k6=l,3 
do 3000 k9=l, 3 

w innn > do 1000 *10=1 

do 1000 kll=l,3 

do 1000 k!2=l, 3 

do 1000 kl3=l, 3 

do 1000 kl4=l,3 

do 1000 k!5=l, 3 

do 1000 kl6=l, 3 

do 1000 kl7=i r 3 

do 1000 kl8=l, 3 
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do 1000 kl9=l,3 
do 1000 k20=l,3 



nbase ( 1 ) = 


V, 1 

Kl 


nbase ( 2 ) — 


KZ 


nbase ( 3 ) = 


ir i 
Kj 


nbase ( 4 ) = 


K H 


nbase ( 5 ) — 




nbase ( 6 ) = 


K D 


nbase ( / ) — 


K 1 


nbase ( 8 ) = 


U Q 

Ko 


nbase ( 9 ) - 




nbase ( i u j 


— n n 


nbase (11) 


= kll 


nbase ( 12 ) 


= kl2 


nbase (13) 


= k!3 


nbase { 14 ) 


= kl4 


nbase ( 15 ) 


= kl5 


nbase (16) 


= kl6 


nbase ( 17 ) 


=kl7 


nbase (18) 


-kl8 


nbase (19) 


= kl9 


nbase (20) 


= k20 



1 

2 
3 



1200 
c 



1250 



do 1250 nn=l, j j 
n=0 

do 1200 j=l,nsub 

if (mset (nn, j ) . eq . 
mset ( nn, j ) . eq . 
mset (nn, j ) . eq . 
mset (nn, j ) . eq , 
n=n+ 1 
endi f 
continue 



if (n.lt . ndif f ) 
goto 1000 
endi f 

continue 



then 



1 .and. nbase (j ). ne . 1 .or. 

2 .and. nbase ( j ) . ne . 2 .or. 

3 .and. nbase ( j ) . ne . 3 .or. 

4 .and. nbase ( j ) . ne . 4 ) then 



1100 



j j=j j + l 

write ( * , 130) (nbase (i ) , i = l, nsub) f j j 
do 1100 i=l,nsub 

mset ( j j , i ) =nbase ( i ) 
cont inue 



1000 
c 

130 



120 
c 



cont a nue 
write ( * , * ) 

format (10x,20(lx,il),5x,i5) 
write ( * , * ) 
write!*, 120) jj 

format (ix, 'Number of words=',i5) 



end 
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SEQUENCE LISTING 

( I : GENERAL I N FORMATION : 

(i> APPLICANT: David W. Martin, Jr. 

^Iclt^^^ Measurement of Gene Expression profiles in 

(lii) NUMBER OF SEQUENCES: 7 
(iv) CORRESPONDENCE ADDRESS * 

i'e! ^f^^f'Kj'jeSte^Kir' ^ Th «*P«"ic S . inc. 

(C) CITY: Hayward 

(D) STATE: California 
(EJ COUNTRY: USA 

(F) ZIP: 94545 

(v) COMPUTER READABLE FORM * 

(A) MEDIUM TYPE: 3.5 inch disk-rtP 

(B) COMPUTER: IBM compatible 

CJ OPERATING SYSTEM: Windows 3 1 
(D) SOFTWARE: Microsoft Word 5 i 

fc (vi) CURRENT APPLICATION DATA* 

(A) APPLICATION NUMBER- 

(B) FILING DATE: 

(C) CLASSIFICATION: 

{vii) PRIOR APPLICATION DATA * 

(A) APPLICATION NUMBER : PCT/US96/095 1 ^ 
(Bi FILING DATE: 06-JUN-96 3 

(vii; PRIOR APPLICATION DATA • 

(A) APPLICATION NUMBER ; PCT/US95/127Q1 

(B) FILING DATE: 12-OCT-95 

(viii) ATTORNEY /AGENT INFORMATION* 

(A) NAME: Stephen C. Macevicz 

(B) REGISTRATION NUMBER • 30 285 

(C) REFERENCE / DOCKET NUMBER : 813wo 

(ix) TELECOMMUNICATION INFORMATION * 

(A) TELEPHONE: (510) 670-936^ 
CBJ TELEFAX: (510) 670-9302 " 



(2J INFORMATION FOR SEQ ID NO: 1 

(i) SEQUENCE CHARACTERISTICS - 

(A) LENGTH: 11 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 
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(xii SEQUENCE DESCRIPTION: SEQ ID NO: 1: 



CTAGTCGACC A 



(2) INFORMATION FOR SEQ ID NO: 2: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 nucleotides 
t'B) TYPE: nucleic acid 
(C) STRANDEDNESS: single 
{ D ) TOPOLOGY : 1 i near 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 



NRRGATCYNN N 



(2) INFORMATION FOR SEQ ID NO: 3: 



(i) SEQUENCE CHARACTERISTICS: 

{A) LENGTH: 38 nucleotides 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
CD) TOPOLOGY: linear 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 



GAGGATGCCT TTATGGATCC ACTCGAGATC CCAATCCA 36 



(2) INFORMATION FOR SEQ ID NO: 4: 



fi; SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 nucleotides 
(3) TYFE: nucleic acid 
{C} STRANDEDNESS: double 
(D) TOPOLOGY: linear 



(xii SEQUENCE DESCRIPTION: SEQ ID NO: 4: 



AGTGGCTGGG CATCGGACCG 2 0 



;2i INFORMATION FOR SEQ ID NO: S: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 20 nucleotides 
IB) TYPE: nucleic acid 
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<C) STRANDEDNESS: double 
(D) TOPOLOGY: linear 

fxi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 
GGGGCCCAGT CAGCGTCGAT 

12) INFORMATION FOR SEQ 10 NO: 6: 

(i) SEQUENCE CHARACTERISTICS • 

,B ' riPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

<*i) SEQUENCE DESCRIPTION: SEQ ID NO: 6: 
ATCGACGCTG ACTGGGCCCC 

(2) INFORMATION FOR SEQ ID NO: 7: 

Ji) SEQUENCE CHARACTERISTICS • 

' TYPE: nucleic acid 
iC) STRANDEDNESS: double 
(D) TOPOLOGY : linear 

txii SEQUENCE DESCRIPTION: SEQ ID NO: 7 : 

AAAAGGAGGA GGCCTTGATA GAGAGGACCT GTTTAAACGG AT^TTCr 

TCTTCC?C' T>r? CC " * w 
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I claim: 

1 . A method of determining the toxicity of a compound, the method comprising 
the steps of: 

administering the compound to a test organism; 

extracting a population of mRNA molecules from each of one or more tissues 
of the test organism; 

forming a separate population of cDNA molecules from each population of 
mRNA molecules from the one or more tissues such that each cDNA molecule of a 
separate population has an oligonucleotide tag attached, the oligonucleotide tags 
being selected from the same minimally cross-hybridizing set; 

separately sampling each population of cDNA molecules such that 
substantially all different cDNA molecules within a separate population have different 
oligonucleotide tags attached; 

sorting the cDNA molecules of each separate population by specifically 
hybridizing the oligonucleotide tags with their respective complements, the respective 
complements being attached as uniform populations of substantially identical 
complements in spatially discrete regions on one or more solid phase supports; 

determining the nucleotide sequence of a portion of each of the sorted cDNA 
molecules of each separate population to form a frequency distribution of expressed 
genes for each of the one or more tissues; and 

correlating the frequency distribution of expressed genes in each of the one or 
more tissues with the toxicity of the compound. 

25 2. The method of claim 1 wherein said oligonucleotide tag and said complement 
of said oligonucleotide tag are single stranded. 

3. The method of claim 2 wherein said oligonucleotide tag consists of a plurality 
of subunits, each subunit consisting of an oligonucleotide of 3 to 9 nucleotides in 

30 length and each subunit being selected from the same minimally cross-hybridizing set. 

4. The method of claim 3 wherein said one or more solid phase supports are 
microparticles and wherein said step of sorting said cDNA molecules onto the 
microparticles produces a subpopulation of loaded microparticles and a subpopulation 

35 of unloaded microparticles. 

5. The method of claim 4 further including a step of separating said loaded 
microparticles from said unloaded microparticles. 
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6. The method of claim 5 further including a step of repeating .aid steps of 
*mphng, sortmg. and separating mtil . nun)ber of Mid ^ ^ ^ ° 

accumulated is at least 10.000. Fancies is 

5 

leas, ,«T;r 0d ° f C ' aim 6 WhCrei " ^ -'Nicies is a, 

■ 0 L, sooner 0 " ° f Cl3im 7 Wherei ° nUmbCr ° f — - a. 

9. The method of claim 5 further including a step of repeating said steps of 
sampl.ng, sortmg, and separating „„„, . nU mber of said loaded mLp art X i 
accumulated ,s sufficient to estimate the relative abundance of a cDNA molecule 

M/. confidence hmtt no larger than 0.1% of said population. 

1 0. The method of claim 4 wherein t~e.. 

,„,,„„ hCrein Sa,d test or ganism is a mammalian tissue 



15 



20 



0 



* culture 



hI P a,oc^s. me,h0d ° f daim ' ° Wherei " ^ " *~ «*« <~ 

12. The method of claim 4 wherein saiH t»c» ~- 

th( >„ wnerein said test organism is an animal selected from 

•he group cons.st.ng of rats. m.ce. hams.er, guinea pig, rabbi,, cats, dog, pigs. ™ 

1 3. The method of claim 12 wherein said one or more tissues are selected from the 
group con-™, of liver, kidney, brain, cardiovascular, thyroid, spleen aTenlle 
•ntesnne sma.l intestine, pancrease „ r ,„ary bladder, stomach, ova.. ,e tes Z d 
mesenteric lymph nodes. 



M. A method of identifying genes which are differentially expressed in a selected 
u-« of a - animal after treatment with a compound, the metho'd 



administering the compound to a test animal; 
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extracting a population of mRNA molecules from the selected tissue of the 
test animal; 

forming a population of cDNA molecules from the population of mRNA 
molecules such that each cDNA molecule has an oligonucleotide tag attached, the 
5 oligonucleotide tags being selected from the same minimally cross-hybridizing set; 

sampling the population of cDNA molecules such that substantially all 
different cDNA molecules have different oligonucleotide tags attached; 

sorting the cDNA molecules by specifically hybridizing the oligonucleotide 
tags with their respective complements, the respective complements being attached as 
1 0 uniform populations of substantially identical complements in spatially discrete 
regions on one or more solid phase supports; 

determining the nucleotide sequence of a portion of each of the sorted cDNA 
molecules to form a frequency distribution of expressed genes; and 

identifying genes expressed in response to administering the compound by 
1 5 comparing the frequencing distribution of expressed genes of the selected tissue of the 
test animal with a frequency distribution of expressed genes of the selected tissue of a 
control animal. 

- 15. The method of claim 14 wherein said oligonucleotide tag and said 
20 complement of said oligonucleotide tag are single stranded. 

16. The method of claim 15 wherein said oligonucleotide tag consists of a 
plurality of subunits, each subunit consisting of an oligonucleotide of 3 to 9 
nucleotides in length and each subunit being selected from the same minimally cross- 

25 hybridizing set. 

1 7. The method of claim 16 wherein said one or more solid phase supports are 
microparticles and wherein said step of sorting said cDNA molecules onto the 
microparticles produces a subpopulation of loaded microparticles and a subpopulation 

30 of unloaded microparticles. 

18. The method of claim 17 further including a step of separating said loaded 
microparticles from said unloaded microparticles. 

35 19. The method of claim 1 8 further including a step of repeating said steps of 
sampling, sorting, and separating until a number of said loaded microparticles is 
accumulated is at least 1 0,000. 
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leas, 1 00 T ";00 me,hOdOrClaim '^"^of.oadedm.cro Panic.es isat 

L 5oIl meth ° d ° f di>im 20 WhCTein is a, 

san.pl.ng, sorting, ^ ^ a numbe P mi lXS ■ 

pnesen, , sa,d p„ pulatlon „ , frequency wj|h . n ^ 

95/.co„nd e , Ke | lm i tno i a ^ erlhan010/oofsa . dpopu|ai . on o5/.w„ha 
f„ ' , ^ r""" 1 ° f Cla ' m ' 7 Wherei " Mimal <* from the group 

24. The method of claim 23 wherein said selected .issue is selected from the 
group cons.st.ng of liver, kidney, brain, cardiovascular thyroid spleen 7n 
intestine smal, intesUne. pancrease u„„a^ bidder, storn^ v^ st. I " 
- mesenteric lymph nodes. ' 

the .^IT ttChniqUe ° f ParalIel "H— -1"™*, .o determine 

the ,o*,c«y of a compound in a test organism, the use comprising the steps of 
adm.rustermg the compound to a test organism 

of the ,eT C ' in8 " P0PUla,i ° n ° f mRNA m ° leCUleS «» «-» - — tissues 

or i:::,' :r sm md > ° f ^ ■»■«*. *>. «* of :„i 

m„.e , de '7 in r *' nUC ' e °" de SeqUenK ° f a P 0 "' 0 " of the cDNA 

modules of each separate population using massively paralle. sig„a, U re sequencing 
to form a frequency distribution of expressed genes for each of rtj SK,uencl "e 
0 .issues; and ™ i enes •<>' "ch of the one or more 



correlating .he frequency distribution of expressed genes in each of ,h, „ 
more tissues with the toxicity of .he compound. ^enes ,„ each of the one or 

26. The use of claim 25 wherein said ,es, organism is a mammahan tissue culture. 

z^:r of cUim 26 wherein ^ — v ™ *- «*« -p— 
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28. The use of claim 25 wherein said test organism is an animal selected from the 
group consisting of rats, mice, hamsters, guinea pigs, rabbits, cats, dogs, pigs, and 
monkeys. 

5 29. The use of claim 28 wherein said one or more tissues are selected from the 
group consisting of liver, kidney, brain, cardiovascular, thyroid, spleen, adrenal, large 
intestine, small intestine, pancrease urinary bladder, stomach, ovary, testes, and 
mesenteric lymph nodes. 

10 30. A use of the technique of massively parallel signature sequencing to identify 
genes which are differentially expressed in a test organism after treatment with a 
compound and which are correlated with toxicity of the compound, the use 
comprising the steps of: 

administering the compound to the test organism; 
1 5 extracting a population of mRNA molecules from a selected tissue of the test 

organism and forming a population of cDNA molecules; 

determining the nucleotide sequence of a portion of each of the cDNA 
molecules using massively parallel signature sequencing to form a frequency 
* distribution of expressed genes; 
20 identifying genes expressed in response to administering the compound by 

comparing the frequencing distribution of expressed genes of the selected tissue of the 
test organism with a frequency distribution of expressed genes of the selected tissue 
of a control organism; and 

determining whether the genes expressed in response to administering the 
25 compound are correlated with toxicity of the compound in the test organism. 



RMsnoom <wo <*7i:w77Ai i > 
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